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ABSTRACT 

Several models for the hard X-ray Background (XRB) suggest that it is due to the 
emission from heavily obscured AGN. Recent studies have revealed the presence of a 
new population of hard X-ray sources which must contribute significantly to the XRB. 
However, whether the majority of these sources are obscured AGN or some other type 
of object at present remains unclear. Here, we further examine the possibility that 
a significant fraction of the XRB comes from a population of galaxies undergoing 
advection-dominated accretion in the high-m regime and thus produce intrinsically 
hard spectra. When the accretion rate is close to m cr ;t, above which an advection- 
dominated accretion flow (ADAF) no longer exists, the major contribution to X-ray 
emission is due to inverse Compton scattering of the soft seed photons produced by 
cyclo-synchrotron radiation. In this regime, the resulting ADAF spectra are relatively 
hard with a fairly constant X-ray spectral index a ~ 0.2 — 0.4 and a spectral cut-off 
at ~ 200 keV. We show that the integrated emission from such sources can provide 
a good fit to the hard (> 2keV) X-ray background, provided that the spectrum is 
dominated by the contribution from objects located at redshifts z ~ 2 — 3. The model 
requires most of the contribution to the XRB to be due to objects accreting at m cr it. 

Key words: galaxies: active - galaxies: nuclei; accretion, accretion disks - X-rays: 
general - galaxies 
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1 INTRODUCTION 

The origin of the X-ray background (XRB) is still largely 
an unsolved problem. Recent progress indicates that the 
XRB is probably due to contributions from different types of 
discrete sources. Although Active galactic Nuclei (Seyferts 
and Quasars, hereafter AGN) provide a large fraction of the 
soft XRB below ~ 2 keV (Shanks et al. 1991; Hasinger et 
al. 1998; Schmidt et al. 1998) their spectra are too steep to 
explain the hard (> 3keV) XRB. Also, as a larger sample of 
QSOs became available, detailed studies of the QSO X-ray 
luminosity function (Boyle et al. 1994) and the source num- 
ber count distribution (Georgantopoulos et al. 1996) have 
shown that QSOs are unlikely to form more than 50 per 
cent of the XRB, even at 1 keV. Models based on the uni- 
fied Seyfert scheme, in which a large fraction of the emis- 
sion from AGN is absorbed by obscuring matter (Madau et 
al. 1994; Celotti et al. 1995; Comastri et al. 1995) seem to 
account for the amplitude and the general shape of the hard 
XRB. However, in these models it is not clear whether cos- 
mological evolution can ultimately smooth out the predicted 
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absorption edges and emission lines (e.g. Gilli et al. 1998) in 
the 2 — lOkeV spectral range, produced by iron and other 
metals (Matt & Fabian, 1994), to a level which is consistent 
with the XRB spectrum statistical errors. In fact, recent 
ASCA observations (e.g. Gendreau et al. 1995) showed that 
the spectrum of XRB is typically very smooth and well ap- 
proximated by a power law of spectral index a — 0.4. Thus, 
to explain the origin of XRB we might have to postulate 
the existence of a new population of hard X-ray emitting 



Recent deep ROSAT surveys have begun to resolve such 
sources at the faintest soft-X-ray flux limits (Vikhlinin et 
al. 1995; Boyle et al. 1995; Almaini et al. 1996; Hasinger et al. 
1998; McHardy et al. 1998). Many of these sources appear 
to be associated with faint active galaxies with implied X- 
ray luminosities 100 times higher than normal field galaxies. 
Furthermore, a promising population of sources with hard 
X-ray spectra has been discovered by ASCA and BeppoSAX 
in the 2-10 keV band (Uedaet al. 1998; Giommi et al. 1998). 
These findings indicate the presence of a population which 
could significantly contribute to the XRB (resolving ~ 30 
per cent of the 2 — lOkeV XRB). Although it seems plau- 
sible that a fraction of the observed hard-spectrum sources 
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might simply be heavily obscured AGN, in the majority of 
cases the true nature of the emission remains unclear. In ad- 
dition, most models predict that the bulk of the hard XRB 
originates at redshifts z ~ 2 which are currently largely un- 
explored in X-rays. 

Continuing on the lines of Di Matteo & Fabian (1997; 
hereafter DF), in this paper we further examine the possi- 
bility that the sources contributing to the XRB have intrin- 
sically hard X-ray spectra, such as are naturally produced 
in the context of advection-dominated accretion models. 

In an advection-dominated accretion flow (ADAF) most 
of the energy released by viscous dissipation is stored within 
the gas and advected inward with the accreted plasma; only 
a small fraction is radiated away (see Narayan, Mahade- 
van & Quataret 1998 for a recent review). Recent work 
on ADAFs has concentrated on low-M solutions (Ichimaru 
1977; Rees et al. 1982; Narayan & Yi 1994, 1995; Abramow- 
icz et al. 1995) which exist when the accretion rate is lower 
then a critical value, m cr it ~ 1.3a 2 , in Eddington units. This 
optically-thin branch is based on two critical assumptions: 
(1) viscous energy is deposited into the ions and only a small 
fraction of energy goes directly into the electrons; and (2) the 
energy transfer between the ions and electrons occurs only 
via Coulomb collisions. Because Coulomb transfer in the the 
low-density ADAF is inefficient, as a result of these assump- 
tions, the accreting gas takes the form of a two-temperature 
plasma, in which the ions are at nearly virial temperature 
T\ = 10 12 /rK and the electron temperature saturates at 
around 10 9 - 10 10 K. 

The spectrum from an optically-thin ADAF is very dif- 
ferent from that of a standard thin disk. The high electron 
and ion temperatures and the presence of an equipartition 
magnetic field allow a variety of radiation processes to con- 
tribute to the emitted spectrum from radio to gamma rays. 
The radio to X-ray emission is produced by the electrons 
whereas the gamma-ray emission is produced by the ions 
via pion production. Here, we restrict the discussion to elec- 
tron emission which comprises synchrotron, inverse Comp- 
ton, and bremsstrahlung processes and in particular to those 
processes which contribute to the X-ray emission. 

In this Letter, we consider whether a significant fraction 
of the hard XRB can be due to the integrated emission from 
a population of galaxies undergoing advection-dominated 
accretion with rh ~ m cr it- When the accretion rate rh is 
substantially lower than rh C rit the X-ray emission from an 
ADAF is dominated by bremsstrahlung radiation. As the 
accretion rate rh approaches m cr it, Comptonized emission 
contributes more significantly to the the high energy spec- 
trum. In a previous paper DF have analyzed the contribution 
to the XRB from ADAFs in the bremsstrahlung-dominated 
regime. Here we show that for sources undergoing accretion 
close to rh ~ m cr it, thermal Comptonization is highly satu- 
rated, and produces hard spectra which can easily account 
for the XRB (this condition then resembles that proposed 
on ad hoc basis by Zdziarski, 1988). 

In §2 of this paper, we compute the detailed spectrum of 
an ADAF in the high-m regime using the model described 
by Narayan, Barret & McClintock (1997) and Esin, McClin- 
tock & Narayan (1997). In §3, we use an appropriate cosmo- 
logical model to derive the contribution to the XRB intensity 
from these sources and compare the resulting spectrum with 
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Figure 1. ADAFs accreting rh < rh cr i t . The seven solid lines 
show the spectra for rh = 0.005, 0.02, 0.04,0.06, 0.07, 0.08,0.089(= 
™crit)i l°g( r tr) = 2.7, a = 0.25, fi = 0.5. The different spectral 
components arc shown for illustration only for the case of rh = 
0.07: the synchrotron and bremsstrahlung emission - dotted line, 
the blackbody emission from the outer thin-disk and the reflection 
component - dash-dotted line, the sum of the Comptonization of 
the soft disk blackbody and synchrotron photons - dashed line. 
The central black hole mass is m = 10 7 . 

the data. We discuss possible implications of our model and 
summarize our results in §4. 



2 MODELING THE X-RAY EMISSION FROM 
AN ADAF 

The hot, two-temperature, advection-dominated solution 
on which our model is based is known to exist only for ac- 
cretion rates below the critical value m cr it (e.g. Narayan & 
Yi 1995; Abramowicz et al. 1995). For rh > rh cr i t Coulomb 
coupling becomes very efficient, allowing the protons to cool 
rapidly and causing the hot flow to collapse into a standard 
cooling-dominated thin disk. In this paper we consider a 
regime of accretion close to rh cr i t . 

In our model, the accretion flow is composed of two 
separate zones: an inner ADAF and an outer thin disk. The 
boundary between the two zones lies at the transition radius, 
r tr (hereafter all radii are in Schwarzschild units). Since we 
are interested in the high rh regime, we adopt a rather low 
value of the transition radius in our calculations, setting 
r t i — 10 2 ' 7 . It is worth noting, however, that the ADAF 
spectrum in the X-ray band in practically independent of 
our choice of the transition radius, as long as it lies outside 
~ 100 Schwarzschild radii (see Esin, McClintock & Narayan 
1997). 

We compute the ADAF spectra using the numeri- 
cal procedure described in detail by Esin, McClintock & 
Narayan (1997) and Narayan, Barret & McClintock (1997). 
This procedure combines a fairly sophisticated treatment 
of radiative transfer (Poutanen & Svensson 1996) with 
global dynamical accretion flow solutions (Narayan, Kato, 
& Honma 1997), while solving for the thermal balance of 
the ions and electrons. In our calculations, we take the black 



© 0000 RAS, MNRAS 000, 000-000 



The X-ray Background 3 



hole mass to be m = 10 7 (in solar mass units), and set the 
viscosity parameter to a = 0.25. We also assume that the 
magnetic field in equipartition with thermal pressure, so that 
the ratio between the thermal and total pressure, denoted 
by /3, is equal to 0.5. 

Spectra produced from such an accretion flow, for a se- 
ries of rh in the vicinity of m C rit , are shown in Fig. hi The low 
energy (radio-IR) end of the spectrum is due to self-absorbed 
synchrotron emission from hot electrons in the equiparti- 
tion magnetic field. The quasi-thermal emission from the 
outer thin disk peaks in the NIR/optical band. The soft syn- 
chrotron and blackbody photons inverse Compton scatter off 
the hot thermal electrons in the ADAF and produce the hard 
power law spectrum with a spectral cut-off at ~ 2 — 3fcT c . 
X-rays scattering off the outer thin disk produce a reflection 
bump at a characteristic energy ~ 30 keV. We also include 
bremsstrahlung emission which peaks at hv ~ kT e . 

The relative importance of Comptonization increases 
very rapidly with increasing rh (see Fig.Q). In the high rh 
regime considered here, the characteristic electron scatter- 
ing optical depth r of the ADAF becomes of order unity, 
since t <x rh. This implies that the value of the Compton 
y— parameter is close to unity as well. As a consequence, in- 
verse Compton scattering becomes the most important cool- 
ing mechanism, and the Comptonized tail dominates the X- 
ray spectrum, as shown by the models plotted in Fig. |l[ 

The energy boost experienced by a photon in each scat- 
tering is proportional to the electron temperature. Although, 
with increasing rh more energy is transferred to the electrons 
by Coulomb coupling, inverse Compton cooling also becomes 
more efficient, with the result that T c saturates at around a 
few times 10 9 K. Since the synchrotron radiation (the main 
source of seed photons) is very soft (see Fig.0), the photons 
will be scattered many times before reaching the Wien sat- 
uration limit, giving rise to a smooth power law spectrum. 
The slope (a x ) of the power law is proportional to logr, 
so that the X-ray spectrum becomes flatter (in L v ) as rh 
increases. 

The hard X-ray spectrum of an ADAF closely resem- 
bles that of the hard X-ray background when the accretion 
rate is in the vicinity of m cr it, but the spectra from sources 
with lower rh become considerably steeper (as illustrated in 
Fig. |l|). However, because the inverse Compton cooling in- 
creases as rh 2 (see Fig. [j]) , eventual contributions to the XRB 
intensity from systems with softer spectra (i.e. smaller rh) 
is not very important, once we assume that a population of 
sources accreting at rh = m cr it exists. Since the sources with 
flat spectra, accreting close to rhcrit, are more luminous, they 
contribute more significantly to the XRB intensity. To make 
this statement more quantitative, the contribution from the 
lower rh regime is negligible, unless the relative fraction of 
the steep-spectrum sources is greater than that for the hard- 
spectrum sources by a factor comparable do the difference 
in their respective hard X-ray luminosities. The latter con- 
dition is satisfied only is we draw our sources from a pop- 
ulation with number density per logarythmic interval in rh 
described by N(rh) oc m _! , where i > 2. 

In the next section we will consider the integrated 
emission from sources producing hard Comptonized spec- 
tra in the high rh regime consistent with a two temperature 
ADAF, and fold it with the appropriate cosmological model 
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Figure 2. The spectrum of the X-ray background (Gruber's em- 
pirical formula indicated by crosses) compared to our best fit 
model spectrum calculated with advection-dominated flows ac- 
creting at rh ~ rh CI .; t and cosmological evolution parameters fixed 
at z max = 3, z = 2, and p = 0. 



to constrain their potential contribution to the unresolved 
3- 100 keV XRB. 



3 CONTRIBUTION TO THE XRB 
3.1 Cosmological Model 

In order to obtain a model for the XRB, we consider the 
contribution from many unresolved sources accreting at 
rh ~ m cr it, distributed over a redshift range from the local 
universe (z — zo) to the early universe (z = z maJ[ ). Comoving 
spectral emissivity from a distribution of such objects can be 
written as a product j[E, z] = n(z)L^,(z), where n(z) is the 
comoving number density of X-ray sources, and Le(z) rep- 
resents the specific luminosity of individual sources (see §2). 
In general, both n and Le are functions of redshift, so that 
neither pure luminosity nor pure number density evolution 
are very good assumptions. Instead, we adopt the following 
simple prescription for the redshift evolution of comoving 
emissivity, j[E,z] = jo(E)(l + z) p , where jo(E) is the lo- 
cal (at z — zo) spectral emissivity and p is the evolution 
parameter. 

The total flux received from such objects, calculated in 
the framework of the Robertson- Walker metric, is given by 

w - ±k x £""" (t\Z^ j ° m + z)]dz > (i) 

where go is the deceleration parameter, Hq the Hubble 
constant (we use go = 0.5 and Ho = 65kms -1 Mpc^ 1 ). 
1(E) represents our computed XRB intensity (in units 
of keVs -1 sr _1 cm - keV -1 ). The only free parameter in 
Eqn. (|l|) is the local source number density, n(zo), which we 
determine by comparing 1(E) with the empirical fit to the 
XRB spectrum. 
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Figure 3. Deviation of our models from the analytical approx- 
imation by Gruber (1992). The error bars are a sample of the 
typical one sigma measurement errors at various energies. Our 
best fit model (heavy solid line also shown in Fig. ^|) gives an ex- 
cellent fit to the observed XRB spectrum over most energy bands. 
It is calculated by considering only the sources with the hardest 
spectra (m = m cr i t ) and by setting z max = 3, zo = 2, p = 0. 
Adding a contribution from objects accreting at different fa in 
the range between 0.005 and 77i cr it (as shown in Fig. [j]; dotted 
and dashed lines) steepens the overall spectrum and the fit wors- 
ens at the lower energies (< 5keV). Adding a 20 per contribution 
to the XRB from Seyfert galaxies with a canonical spectral slope 
of a = 0.7 (shown as a dot-dashed line for fa = fa CI i t model and 
as a dotted line for the model with a range of m) further increases 
the excess at soft energies. To illustrate the effects of cosmology, 
we also plot (thinner solid and dashed line) the spectrum com- 
puted for the same model as that shown by the heavy solid line, 
but with different cosmological evolution parameters: z max = 3.5, 
zo = 0, p = 3. 



3.2 Comparison with the data 

An empirical fit to the spectrum above ~ 3keV where 
most of the energy density resides, has been deter- 
mined by Gruber (1992). The energy flux in units of 
keVs -1 sr _1 cm~ 2 keV" 1 is given by: 

7.877£ , - 0M exp(-£/41.13 keV), 

3keV < E < 60keV; 



Iobs(E) 



1652.E 



" 2 -° + 1.754B" ' 7 . 



(2) 



60keV < E < 100 keV. 



Fig. [2] shows a comparison between the data and our best fit 
model for the predicted contribution of the ADAFs to the 
background. Fig. ^| shows the fractional difference between 
our model results and the analytical approximation to the 
XRB in the energy range 3 - 100 keV (Eqn. ||). To judge the 
goodness of fit,we show a few typical measurements errors at 
different energies. The different fits (described in detail be- 
low) , which can be identified from the plot of their residuals 
in Fig. ^, correspond to essentially two different models. 

In one model, corresponding to the dashed and dot- 
ted line in Fig. [], we take into account contributions from 
sources accreting at a range of mass accretion rates (specifi- 
cally, we consider sources whose spectra are shown in Fig. [j] 
contributing in equal proportion) in each successive redshift 



shells. As we argued in §2, the contribution from objects 
with m lower than the values considered here is negligible as 
long as the relative fraction of sources at a given m < m cr it 
does not exceed the number of sources close to rh cr it by a 
factor corresponding to ~ l^fn- ^ ne thinner dashed- 

line corresponds to a model with 2 max = 3.5 and 20 = 0. In 
this case, the normalization of the XRB spectrum requires a 
local comoving number density of n(zo — 0) ~ 10 -4 Mpc -3 
The best fit evolutionary parameter is p = 3 (similar to that 
inferred for AGN), implying that most of the contribution 
to the XRB intensity comes from objects at high redshifts. 
To illustrate this further, we show a model (thicker dashed 
line) in which all of the emission comes from z m&K — 3 and 
zo — 2 and the evolution parameter p is set to be equal 
to 0. In this case, the comoving number density of sources, 
at these redshifts, is consistently higher and of the order 
of n(z = 2) ~ 10~ 2 Mpc~ 3 . The quality of the two fits 
(one with strong evolution and the other with a single pop- 
ulation of sources between 2 = 2 and 3) is very similar. 
For E > 5keV the computed intensity is within 2a of the 
data points, at E < 5keV, it significantly overestimates the 
observed spectrum. However, small amounts of absorption 
(i.e. a column density Ah ^ a few 10 21 cm -2 - much less 
than that required by Seyfert models of the XRB) would 
significantly improve the fit at these low energies, where the 
emission from other extragalactic sources also becomes rel- 
evant. 

If we only consider the emission from the brightest 
sources, those with rh ~ rh cr i t we find that the fit to the low 
energy range of the XRB significantly improves: the com- 
puted intensity now lies within la of the data points (solid 
lines in Figs, g and ^|). The spectra from ADAF sources ac- 
creting at these high m values are very hard and give an ex- 
cellent fit to the XRB when appropriately redshifted. Again, 
we compute the model intensity considering first, a strongly 
evolving population (p — 3) from 2 max = 3.5 to zo = 
(thinner solid line) and second, a short lived, non-evolving 
(p = 0), population of objects from z = 3 to zo = 2 (heav- 
ier solid line; this our best fit model plotted in Figs. [|). The 
normalizations to the XRB yield comoving number densities 
very similar to those obtained above. 

The same set of models described above (for 2 max = 
3.5, Zo = and p — 3), is also calculated considering a 20 
per cent contribution to the 2—10 keV XRB due to AGN 
characterized by a canonical power-law spectrum with en- 
ergy index q x = 0.7. Such a contribution produces a further 
steepening of the spectrum and therefore a greater excess 
at soft energies (3 — 5keV range; dotted and dashed line in 
Fig. ^|). However, the absorption column density required 
to suppress such an excess is still much smaller than that 
required by Seyfert models of the XRB and of the order of 
Ah 10 22 cm" 2 . Note that because of the outstanding is- 
sue that the observed AGN spectrum does not match that 
of the XRB (there is no strong evidence for any hardening 
of the X-ray spectral slope even for AGN identified in deep 
ASCA fields; Boyle et al. 1998), we know that there has to 
be a very significant residual contribution to the 2 — 10 keV 
XRB from sources that we have yet to identify (which we 
suggest might be ADAFs). 

In summary, we have directly shown that a significant 
contribution to the XRB can result from a population of 
galaxies undergoing advection dominated accretion in their 
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nuclei with rh ~ m cr it- In this high-m regime an ADAF 
is luminous enough to satisfy the energy constraints of the 
XRB (in the 2 — 10 keV band the ADAF luminosities are in 
the range of 2 - 8 x 10 42 ergs -1 for 0.07 < rh < m crit , for 
m = 10 7 as shown in Fig. |l|); moreover, an ADAF produces 
a very hard spectrum, through photon-starved Comptoniza- 
tion, which can explain the observed spectrum of XRB. 



4 CONCLUSION AND DISCUSSION 

The main results obtained in this paper can be summarized 
as follows: 

(a) In the framework of the ADAF paradigm it is pos- 
sible to construct a model which reproduces with good ac- 
curacy the XRB spectrum in the range ~ 3 — 100 keV and 
meets all the presently known constraints on number density 
of source population and X-ray spectral characteristics. 

(b) The key feature of the model is the existence of a 
population of ADAFs accreting at m ~ m cr it- In this regime 
thermal Comptonization of cyclo-synchrotron photons gives 
rise to a typically hard X-spectrum (with a x ~ 0.2 — 0.4 ) 
which can naturally fit the 3 — 100 keV XRB. Also, because 
of the small solid angle subtended by the cold material, any 
spectral features due to X-ray irradiation are negligible, and 
the model is consistent with the observed smoothness of the 
XRB spectrum. 

(c) We compute in detail the typical luminosity and 
spectra for such sources. The contribution from the Comp- 
tonized emission from a population of ADAFs is integrated 
over redshift to reproduce the XRB intensity . Good fits are 
obtained for z ma x ~ 3.5 and values for the evolution param- 
eter p ~ 3. Implying that the peak of the XRB production 
is around z ~ 3 — 2 (similar quality of fits are in-fact ob- 
tained by considering a single population of ADAFs between 
Zmax = 3 and zo — 2 and assuming p = 0) The normaliza- 
tion to the XRB intensity yields a comoving number density 
of no ~ 10 -4 Mpc -3 , similar to that of present day Seyfert 
galaxies. 

The requirement that z max ~ 3 suggests that the XRB 
is produced before the main quasar phase, which peaks at 
redshift z ~ 2. If a population of massive black holes was 
assembled very early, by z ~ 3, maybe in dwarf galaxies (e.g. 
M32 contains a black hole of M B h = 2 x 10 6 M q ; Van der 
Marel 1998; Magorrian et al. 1998) it could have undergone 
advection-dominated accretion between z ~ 3 and 2. It is 
possible that during that time the infall of gas into the cen- 
tral black hole was increasing, e.g. due to galaxy mergers or 
to star formation, which allowed the sources to accrete close 
to m C rit- In this scenario, the end of the ADAF phase would 
simply be due to the transition to the active QSO phase 
with rh > m C rit. We further speculate that once the accre- 
tion fuel was depleted, quasars rapidly dropped again into 
the low-m ADAF state. In this subsequent ADAF phase, m 
is much lower and bremsstrahlung emission can provide a 
further energy contribution to the hard XRB (as discussed 
by DF where z ma x was determined to be ~ 2). 

In this paper we have assumed that all systems have 
black holes with a constant mass, m = 10 7 . Since the ADAF 
spectra are rather insensitive to m while X-ray luminosity 
changes linearly with this parameter, adopting a varying m 
prescription will simply give a different value for the comov- 



ing number density of X-ray sources and will not change any 
of the other results in this paper. 

We must note that decreasing electron temperatures in 
the accretion flow by a factor of ~ 2 — 4 would cause the 
spectra to have lower cutoff energies, moving the redshift at 
which the bulk of XRB is produced closer to the QSO death- 
line at z ~ 2. We believe that the present ADAF model is 
robust enough to make such a decrease in the gas tempera- 
ture unlikely. However, some physical mechanism currently 
not included in the model, for example the presence of strong 
winds from the accretion flow (as suggested by Blandford & 
Begelman 1998), might produce such an effect (e.g. Di Mat- 
teo et al. 1998). 

Deep X-ray observations with AXAF and later missions 
should resolve objects at redshifts beyond the peak of the 
XRB production, thereby determining the origin and the 
nature of the hard sources contributing to XRB emission and 
providing clues for understanding the evolution of the fueling 
of massive black holes. If ADAFs do contribute significantly, 
the absence of a strong UV bump in the spectrum means 
that the objects will not have a normal broad-line region. 
They will not resemble a standard, broad-line AGN in the 
optical/UV band. 

The main limitation of the high-m ADAF model (as it 
has become apparent in §3), is that it requires most of the 
objects to be accreting close to m cr it- Although this might 
seem a peculiar feature it is worth noticing that models for 
the galactic black hole candidate Cyg X-l (Esin et al. 1998), 
can also be explained by an ADAF which spends most of its 
time accreting at m cr it ; this might indicate that mechanisms 
might exist which lock the accretion rate in some systems 
at values close to their critical rate. 
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